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COMPUTING WITH QUADRATIC FORMS 
OVER NUMBER FIELDS 

PRZEMYSLAW KOPROWSKI AND ALFRED CZOGALA 



Abstract. This paper presents fundamental algorithms for computational 
theory of quadratic forms over number fields. In the first part of the paper, we 
present algorithms for checking if a given non-degenerate quadratic form over 
a fixed number field is either isotropic (respectively locally isotropic) or hyper- 
bolic (respectively locally hyperbolic). Next we give a method of computing 
the dimension of the anisotropic part of a quadratic forms. Finally we present 
algorithms computing the level of a number field and verifying whether two 
number fields have isomorphic Witt rings (i.e. are Witt equivalent). 



1. Introduction 



L^ The algebraic theory of quadratic forms is an important branch of mathematics, 

*~^ which has already achieved its maturity. Yet still, the computational side of this 

theory is seriously under-developed. Little has been done so far to develop algo- 

i -^ rithms for fundamental problems on quadratic forms. There are a few papers that 

treat computational aspects of the theory but they mostly concentrate on forms 
over the rationals (see e.g. HHj)- While the algebraic theory of quadratic forms 
over number fields (i.e. finite extensions of Q) are very like the theory over the 
rationals, the computational approach seems to be rudimentary here. The aim of 
this article is to partially fill this gap, as well as provoke farther discussion and 

qq future research. 

^^ This paper is organized as follows: in Section [2] we present an algorithm (see 

Algorithm fjh checking if a given form (over a fixed number field K) is isotropic. 

**~\ This algorithm uses sub-procedures (Algorithms^andpj) deciding whether the form 

is isotropic at a non-archimcdcan prime of K (respectively odd or even) . These two 
algorithm may be of an independent interest to the reader. Next, in Section [3] we 
show Algorithm [8] determining if a quadratic form is hyperbolic. It is known that 
• • any non-degenerate form can be uniquely decomposed into an orthogonal sum of 

its anisotropic part and a hyperbolic form (one of these two parts may of course be 
void if the form in question is either anisotropic or hyperbolic itself) . Half of the 
dimension of the hyperbolic part is an important invariant of the quadratic form, 
known as the Witt index. In Section HI we shows a procedure effectively computing 
the Witt index of a form. 

In Sections [5] and [6] we go a step further and develop algorithm computing in- 
variants of the ground fields, that play important roles in the algebraic theory of 



quadratic forms. Algorithm 12 computes the level s{K) of a number field K, which 
is the length of the shortest representation of — 1 as a sum of squares. 

Recall that the set WK of similarity classes of non-degenerate symmetric bilinear 
forms over a given base field K is a ring with operations induced by the orthogonal 
sum and the tensor product. It is called the Witt ring of the field K. Because 
a bilinear form defines an orthogonal geometry on the vector space on which it is 
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2 PRZEMYSLAW KOPROWSKI AND ALFRED CZOGALA 

defined, thus the Witt ring can be viewed as an algebraic structure encoding infor- 
mation on all possible orthogonal geometries over a given base field. Two fields are 
said to be Witt equivalent, if their Witt rings are isomorphic. In [12| K. Szymiczek 
described the set of global field invariants fully determining its Witt equivalence 
class. In Section [6] we present Algorithm [14] computing all these invariants. In 
particular the algorithm may be used to verify whether two number fields are Witt 
equivalent. 

The authors implemented all the algorithms presented in this paper in a com- 
puter algebra system Sage |llj . Using this implementation, we were able to find 
representatives of Witt classes of number fields of low degrees. The results are pre- 
sented in Tables [T] [3j Moreover, we used our implementation to empirically analyze 
the distribution of different Witt classes among number fields of a fixed degree. The 
quantitative results are gathered in Tables [4}{6j 

In Sections [2] through HI K = Q(#) is always a fixed number field specified by 
the minimal polynomial of $ over Q and Ok is the integral closure of Z in K. In 
most of the algorithms presented in this paper, we need to determine whether a 
certain elements of a number field are squares. It is obvious that o is a square if 
and only if x 2 — a factors in AT[x]. The last condition can be verified using |15[ 
§5.4]. In order to simplify frequent references we write the algorithm explicitly. 

Algorithm 1. Given a element a of a number field K, this algorithm returns true 
if and only if a is a square, otherwise it returns false. 

(1) Use [13 §5.4] to check if x 2 — a factors over K, if so return true, if not then 
return false. 

Similarly we often need to check if an element is a square in the completion K v 
of the number field K at some finite prime p. We present two separate procedures 
respectively for odd and even primes. 

Algorithm 2. Given an element a of a number field K and an odd prime p. This 
algorithm returns true if and only if a is a square in Kp. Otherwise it returns false. 

(1) Compute the valuation ord p a. If it is odd then return false and quit. 

(2) Use [TJ Algorithm 4.8.17] to find a uniformizer 7r of p. 

(3) Check if the class of a ■ it~ mdv a is a square in the residue field K/p. Return 
true if it is a square and false otherwise. 

The proof of correctness of this algorithm is straightforward. Listing [2] contains 
the pseudo-code. The case of even primes is a bit more tricky. 

Algorithm 3. Let d be an even prime of K and a £ Ok be a non-zero element. 
This algorithm returns true if and only if a is a square in the completion Kj,, 
otherwise it returns false. 

(1) Use Algorithm n] to check if a is a square in K, if so return true and quit. 

(2) Take L = K{^/a). Use H Algorithm 2.4.13] to decompose in O l . If 5 
splits into two primes in Ol, then return true, otherwise return false. 

Proof of correctness. It is clear that if x is already a square in K, then it is a 
square in K?, as well. This justifies step (1). Suppose that it is not a square in K, 
take L = K{yfa) and let D be a prime of L dominating t). Assume that splits 
in L, it follows that both: the relative ramification index e( s /o) and the relative 
inertia degree /(^/o) are equal 1 and so is the local degree (L® : if ). Consequently, 
s/a £ L c Lj> = Kjj, as desired. 

Conversely, assume that a is a square in K a . Let S) be any prime of L laying 
over 0. Then L® = K (y/a) by [9J Theorem 5.5], hence Ljj = if and so e^/o) = 
/( s /°) — 1- This shows that, for every prime S of L over 5 we have (Ljj : K ) = 1. 
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But (L: K) = 2 is the sum of the local degrees (L® : A ) for all the primes 2) 
dominating 0. It follows that there are exactly two such primes. □ 

Listing [3] contains the pseudo-code of the algorithm. 

Observation 1.1. There exists an algorithm which, for a non- degenerate quadratic 
form presented as a symmetric matrix, finds a list {ai, . . . , a^} with entries in Ok 
and such that the form (01, . . . , ad) is isometric to <p> over K . 

Proof. The Gram-Schimdt orthogonalization algorithm is well known. Once the 
form is diagonalized, with entries expressed using the power basis, it is straight- 
forward to clear denominators using the fact that a /b and a ■ b belong to the same 
square class. □ 

From now on we will always assume that all quadratic forms are non-degenerate 
and, as inputs to Algorithms Hl-fTTI they are given as lists of entries from Ok- 



2. ISOTROPY OF A QUADRATIC FORM 

In this section, we present an algorithm that checks if a given form <p over a 
number field K is isotropic or not. The organization of this section reflects the 
general idea of solving the problem locally. Hence, Algorithm [5] checks if the form 
is isotropic at an odd prime, farther reducing the task to the residue field (see Algo- 
rithniffl). Next, Algorithm ^verifies whether the form is isotropic at an even prime. 
Finally, Algorithm [7] checks if the form is globally isotropic, using the two above- 
mentioned algorithms as sub-procedures. Recall (see e.g. [T2J Definition 15.2.1]) 
that the discriminant of a quadratic form </? is defined by the formula 

disc( /? :=(-l) d(d " 1, / 2 det(^, 

where d = dim (p. 

Algorithm 4. Let p be an odd prime of K and ip = (01, . . . , a<j) be a quadratic 
form with all its entries being p-adic units. This algorithm returns true if and only 
if the residual form tp ® K/p is isotropic, otherwise it returns false. 

(1) If dimip = 1, return false. 

(2) If dinup = 2, return true when disctp is a square in K/p, otherwise return 
false. 

(3) If dimip > 2, return true. 

The correctness of the above algorithm follows immediately from [3 Theo- 
rem 1.3.2]. Listing B] contains the pseudo-code of the algorithm. 

Algorithm 5. Let p be an odd prime of a number field K. Given a quadratic form 
ip, this algorithm returns true if ip £§) K 9 is isotropic and false otherwise. 

(1) If dump = 1, return false. 

(2) If dim tp > 5, return true. 

(3) Let {ai, . . . , ad] be the list of coefficients of ip, all a, € Ok- Partition this 
list into two sublists depending on the parity of the p-adic valuation: 

Lpo := {a, ■ TT- m ' d " Ql J ordp a. t = (mod 2)}, 

<pi := {a 7 . ■ 7r-° rd » Q * J ordp a t = 1 (mod 2)}. 

Here n is a uniformizer of p computed using jTJ Algorithm 4.8.17]. 

(4) Use Algorithm p] to verify whether any of <po, <px is isotropic over A"/p. 
Return true if Algorithm [4] returned true at least once, otherwise return 
false. 
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The correctness of the algorithm follows from [7J Proposition VI. 1.9]. The 
pseudo-code of the algorithm is presented in Listing [5j Now, we consider even 
primes. Recall (see e.g. [7J Definition V.3.17]) that the Hasse invariant of a qua- 
dratic form if = (oi, . . . , a<j) at a prime p is: 

(1) s p (<p):= Yl ( a i> a i)p> 

l<i<j<d 

where (cii,aj)p denotes the p-adic Hilbert symbol. Recently J. Voight in [14] pre- 
sented an algorithm for computing the Hilbert symbol in a completion of a number 
field. We use it to verify whether a quadratic form is isotropic over a dyadic com- 
pletion of K. 

Algorithm 6. Let be an even prime of K and <p be a quadratic form over K. 
This algorithm returns true if and only if p ® Kx, is isotropic, otherwise it returns 
false. 

(1) If dim ip < 1, then return false and quit. 

(2) If dimip = 2, then use Algorithm [3] to check whether disct^ is a square 
in Kjj. If so, then return true and quit, otherwise return false and quit. 

(3) If dimtp = 3, then proceed as follows: 

(a) Use (14,, Algorithm 6.6] to compute the Hilbert symbol (—1, — det (<£>)). 

(b) Use Eq. pj and [131 Algorithm 6.6] to compute the Hasse invariant 
So(<p) of ip at 0. 

(c) If (— 1, — det(y?)) = sq((p), then return true otherwise return false 

(4) If dimip = 4, then proceed as follows: 

(a) Use Algorithm [3] to check if det p is a square in Kx,. If not, then return 
true and quit. 

(b) If det ip e {K%) 2 , then use Eq. Q and QH Algorithm 6.6] to compute 
the Hasse invariant s$(ip) and the Hilbert symbol (— 1, — 1) . Return 
true if they are equal, return false if they are not. 

(5) If dimtp > 5, then return true. 

Proof of correctness. An unary form is never isotropic and a quintic or higher- 
dimensional form over a dyadic field is always isotropic by the means of Theo- 
rem VI. 2. 12]. This justifies steps [T] and [5] Next, it is well known that a binary form 
is isotropic if and only if its determinant is a minus square, which proves step pi). 
On the other hand, if the form has dimension three, then [7J Proposition V.3.22] 
asserts that it is isotropic if and only if (— 1, — det(ip)) = Sx,(<p). 

This leaves us with quaternary forms. Now, [7J Corollary VI. 2. 15] asserts that 
over a local field there is only one anisotropic form of dimension 4 and its determi- 
nant is a square. Thus, if det ip $. (K^) 2 , then ip Kx, is necessarily isotropic. On 
the other hand, if det ip € (K£) 2 , then [7J Proposition V.3.23] provides us with a 
needed criterion for isotropy □ 

We present the pseudo-code of this algorithm in Listing [61 Now, we are ready 
to present the main algorithm of this section, that checks if a form is isotropic over 
a given number field. 

Algorithm 7. Given a quadratic form ip — (ai , . . . , ad) over K with aj € Ok , this 
algorithm returns true if and only if tp is isotropic and false if it is not. 

(1) If dim tp < 1, then return false and quit. 

(2) If dim^ = 2, then use Algorithm [l] to check if discy is a square in K. If 
so, then return false; if not, return true. 

(3) Use a real-root separation algorithm (e.g. [SJ §8.5.2]) on the generating 
polynomial of K to find all the real embeddings pi, . . . , p r of K. Compute 
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the signature Sj := sgn Pj(p) of image of p under each embedding. If 
\sj\ — dim p for any 1 < j < r, then return false and quit. 

(4) Use 1, §6.2.5] to factor 20 K into prime ideals 20 K = d^ 1 ■ ■ ■ X>%* in O k . 
For each D t use Algorithm [6] to check if p^K^. is isotropic. If the algorithm 
returns false, for at least one t>i, then return false and quit. 

(5) Use !2 ( 2.3.22] to find all odd primes p of K dividing any of the coefficients 
at of p. For each such a prime p call Algorithm [5] If the procedure returns 
false at least once, then return false and quit. 

(6) Return true. 

Proof of correctness. The cases of unary and binary forms are trivial. For forms 
of higher dimension we use the local-global principle 7, Principle VI. 3.1]. The 
form is isotropic over K if it is isotropic over all the completions of K. Now p, 
having dimension at least three, is trivially isotropic at all odd primes that do not 
divide any of the coefficients. These are almost all primes of K. Thus, we are left 
with only finitely many cases to check: finitely many real places treated in step [3j 
finitely many dyadic places covered by step HI and finitely many non-dyadic primes 
considered in step [5] □ 

The pseudo-code of Algorithm [7] is presented in Listing [7j 

Remark 2.1. In order to compute in step [3] the signatures of p under real em- 
beddings of K, one may proceed as follows. Write all the roots i?i, • • ■ , f? r of the 
defining polynomial / in the interval representation as explained in [8j §8.5]. Let 
a.j = a,jQ + aji'd + • • • + aj !n _i# n be a coefficient of p expressed in the power 
basis representation. Take gj :— Ojo -\-dj\X-\- ■ ■ ■ +a,j tn -\X n and use [8J §8.5, Sign 
evaluation] to compute Sji = sgnu/j^i)) for every 1 < i < r. Then Si — J^ Sji. 

3. HYPERBOLICITY OF A QUADRATIC FORM 

In this section we present an algorithm checking another fundamental property 
of a quadratic form, namely whether it is hyperbolic (hence, a zero element in the 
Witt group). The general idea is similar to the one adopted in the previous section. 
Again, we treat the problem locally, separately for odd primes, even primes and 
real embeddings of K. 

Algorithm 8. Given a quadratic form p = {a%, . . . , ad) over K with a, € Ok, this 
algorithm returns true if and only if p is hyperbolic, otherwise it returns false. 

(1) If dim(y9 is odd, then return false and quit. 

(2) Compute the discriminant disc p. Use Algorithm [Tl to check if disc p> is a 
square in K. If it is not, then return false. 

(3) Use a real-root separation algorithm (e.g. [SJ §8.5.2]) on the generating 
polynomial of K to find all the real embeddings p\, . . . ,p r of K . Compute 
the signature Sj :— sgn pj (p) of image of p under each embedding. If Sj =/= 
for any 1 < j < r, then return false and quit. 

(4) Find all odd primes of K dividing any of the coefficients a^ of p, using [51 
2.3.22]. For each such a prime p proceed as follows: 

(a) Let 7r be the uniformizer of p found using [1, Algorithm 4.8.17]. 

(b) Take p x := {o, ■ 7r~ ord p a - \ ord p a % = 1 (mod 2)}. 

(c) If dim pi is odd, then return false and quit. 

(d) If the residue class of disc</?i is not a square in the residue field K/p, 
then return false and quit. 

(5) Factor 2 over Ok using [TJ §6.2.5] to find all even primes of K. For each 
even prime compute the Hasse invariant s$((p) using Eq. (fij) and [141 
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Algorithm 6.6]. If Sg(</?) ^ 1 for at least one prime 0, then return false and 
quit. 
(6) Return true. 

Proof of correctness. It is well known that the discriminant of a hyperbolic form is 
a square and its dimension has to be even. Moreover, the signature of a hyperbolic 
form over any real closed field is zero. Therefore, the form tp over K must have 
signature zero under every real embedding of K, in order to by hyperbolic over K. 

Now, suppose that a form p of an even dimension and with discriminant being 
a square has a zero signature in every ordering of K. Fix an odd prime p and let 
pi := {ai ■ 7r-° rd » a * | ord p a 4 = 1 (mod 2)}. It follows form Q31 Theorem 9.2.3], 
that he second residue homomorphism of <p, with respect to p, is zero if and only if 
the form tpi ® K/p is even-dimensional and its discriminant is a square (in K/p). 

If the second residual homomorphisms with respect to all odd primes of K are 
null, then the Witt class of ip sits in m{WO K ) n I 2 K by Ch. IV, Corollary 
4.5], where ^(WOk) denotes the nilradical of the Witt ring of Ok and IK is 
the fundamental ideal of the Witt ring WK. Clearly one needs to check only 
these primes that divide any of the coefficients of p> as we do in step |4j Now, 
let Di,...,Q g be all the dyadic primes of K, [4j Proposition 3.5] asserts that the 
map <p n- (so 1 (y), . . . , Sj )g _ 1 ((p)) as an isomorphism from %1(WOk) H I 2 K onto 
{±1} 9_1 . Thus, ip is hyperbolic if and only if all its dyadic Hasse invariants vanish. 
This proves the correctness of the algorithm. □ 

In the previous algorithm we did not need a separate sub-procedure for checking 
if a given form is hyperbolic at a given finite prime p of K. Nevertheless, we present 
such an algorithm below. First of all, it nicely mirrors Algorithms [5] and [6] from the 
previous section. More importantly, it is used in the next section. 

Algorithm 9. Let p be a finite prime of a number field K (either even or odd). 
Given a quadratic form ip, this algorithm returns true if the form <p v := ip ® K p is 
hyperbolic and false otherwise. 

(1) If dim 93 is odd, then return false and quit. 

(2) Compute the discriminant disc p> and check if it is a square in the completion 
K p (use either Algorithm [2] or [3l. If it is not a square, then return false 
and quit. 

(3) use Eq. ([I]) and [HI Algorithm 6.6] to compute the Hasse invariant s p (p) 
and the power (— 1, — 1)™ of p-adic Hilbcrt symbol, where 2m = 
dim 95. Return true if they are equal, return false if they are not. 

Proof of correctness. Take a form p of an even dimension. If the discriminant disc ip 
is a square in K v and the Hasse invariant Sp(ip) equals (— 1,— l)p , then <p 

is isometric to the hyperbolic space m(l, —1) by Proposition V.3.25]. □ 



4. Witt index of a quadratic form 

Recall (see e.g. [JJ Chapter i, §4]) that any non-degenerate quadratic form p 
can be uniquely (up to an isometry) decomposed as <p> = ip _L H , where ip is an 
anisotropic form, called the anisotropic part of <p and H is hyperbolic. The number 
of hyperbolic planes constituting H (i.e. half of the dimension of H) is called the 
Witt index of p> and denoted ind(y). In this section we present an algorithm that 
computes the dimension of the anisotropic part of <p. It can be also used to educe 
the Witt index since clearly ind(/3 = 1 /2 ■ (dim 95 — dimi/;). Again, the problem is 



the global one in Algorithm 11 



first solved locally (see Algorithm 10 1 and then the local solution is used to derive 
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Algorithm 10. Given a non-degenerate quadratic form ip over a number field K 
and a finite prime p, this algorithm computes the dimension of the anisotropic part 
of ip p := tp ® K p over the completion K p . 

(1) dim (p is even, proceed as follows: 

(a) Use Algorithm [9] to check if <p p is hyperbolic. If it is, then return 
and quit. 

(b) Check if disc tp is a square in K p using either Algorithm^ (if p is odd) 
or Algorithm^ (when p is even). If so, then return 2 and quit. 

(c) Return 4. 

(2) dim(p is odd, proceed as follows: 

(a) Let n := dimtp and take ip :— ip _L ((— i)" ( " + 1} / 2 . det<p). 

(b) Use Algorithm |9] to check if tp (g) K p is hyperbolic. If it is, then return 
1 and quit. 

(c) Return 3. 

Proof of correctness. First assume that <p is an even-dimensional form, so <p p £ IK p . 
If it is not hyperbolic, then its class in the Witt ring WK p is not zero. Suppose 
that disci/? is a square in K p . It follows that ip p £ I 2 K p . But for a local field there 
is only one non-zero element of I 2 K p , namely the form n p — (1,u,tt,utt), here u 
is a p-adic unit such that K p (y/u) is the unique unramified extension of K p (sec 
[7J Corollary VI. 2. 15]). It follows that the anisotropic part pf <p p has dimension 4. 
Conversely, suppose that disctp is not a square in K p . Therefore ip p £ IK p \ I 2 K p 
and so the anisotropic part of ip p has dimension 2. 

Now assume that the dimension of <p is odd. Hence, the form ip constructed in 



step ( 2a ) is an even dimensional form and its discriminant is a square in K p . Conse- 
quently, the Witt class of i/j p := ip ® K p sits in I 2 K p . If ip p is hyperbolic then i/j p = 
=±1(1,-1), hence <p p _L (1,-1) = (c) _L =±1(1,-1) for c = -(-1)" ( " + 1] / 2 det ip. 
This implies that the anisotropic part of p p is unary. Conversely, suppose that ip p 
is not hyperbolic. As in the first part of the proof, this leads to ipp — Vp m the Witt 
ring WK p . In particular, the Witt classes of tp p and (c, l,u,7r,u7r) are equal. But 
a quintic form over a local field is necessarily isotropic and so it is similar to either 
ternary or unary form. We claim that the unary case is impossible. Indeed, if 

(c, 1, U, 7T, utt) = (x) J- 2(1, —1), 

then square classes of c and x are equal and the Witt cancellation theorem asserts 
that the forms (l,u, rr, wk) and 2(1,-1) are isometric over K p contradicting [71 
Corollary VI. 2. 15]. All in all, (c,1,u,tt,utt) has a ternary anisotropic part and so 
has ipp. □ 

Algorithm 11. Given a non-degenerate quadratic form ip over a number field K, 
this algorithm computes the dimension of the anisotropic part of p. 

(1) Use a real-root separation algorithm (e.g. [SJ §8.5.2]) on the generating 
polynomial of K to find all the real embeddings p\ , . . . , p r of K and compute 
the maximum of the signatures 

N := max \sgn Pl (p)\. 

l<j<r 

(2) If N > 3, then return N and quit. 

(3) Find all primes of K dividing any of the coefficients a* of tp, using jH 2.3.22]. 

(4) Factor 2 over Ok using 1, §6.2.5] to find all even primes of K. 

(5) Let C be the set consisting of all even primes of K and all odd primes 
dividing any of the coefficients of p. 

(6) For every p £ C compute the dimension d p of the anisotropic part of (p®K p 



using Algorithm 10 and let M = max{d p | p £ £}. 
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(7) Return max{M,iV}. 

Proof of correctness. Let tp be the anisotropic part of if. Obviously 

dim'0 = dim tp = | sgnpj(ip)\ (mod 2) 

for any real embedding pj of K. Take N to be the maximum of the signatures 
of ip at all the real places. Now, ip being anisotropic must be anisotropic at some 
place of K, either finite or infinite. Therefore, clearly dim?/; is the maximum of 
the dimensions of the anisotropic parts of the localizations of ip at all the places of 
K . However, if TV > 3, then we do not need to consider the finite primes at all. 
Indeed, if dim-0 > 5, then [Jj Theorem VI. 2. 2] implies that it must be an infinite, 
hence real, place. Therefore in this case dim-0 = N > 5. Similarly, if N — 3 or 
N = 4, then there is a real embedding pj, with sgn pj(ip) = sgnpj(ip) — N and so 
dim-0 > N. On the other hand, it cannot be strictly greater, since otherwise ip 
would have to be anisotropic at some finite place contrary to the already mentioned 
Theorem VI.2. 2]. □ 

5. Level of a number field 

In this section we present an algorithm determining an important invariant of a 
number field, namely its level. Recall that a level of a field K, denoted s(K) is the 
minimal number of terms needed to represent —1 as a sum of squares in K. We set 
s(K) = oo, when —1 cannot be expressed as a sum of squares (i.e. K is formally 
real) . 

Algorithm 12. Given a number field K = Q("i?) specified by its defining polyno- 
mial /, this algorithm computes the level s(K). 

(1) Use Sturm's sequence to check if / has any real roots. If so, then return oo 
and quit. 

(2) Use Algorithm [l| to check if —1 is a square in K. If so, then return 1 and 
quit. 

(3) Use |U §6.2.5] to find the factorization of 2 in Ok in the form of a list C 
consisting of triples (j)j, e^, /,•), where Oj is a prime of K dominating 2 with 
the ramification index ej and the inertia degree fj. 

(4) If for any j, both Cj and fj are odd, then return 4 and quit. 

(5) Return 2. 

Proof of correctness. The standard Sturm's sequence method counts the number 
of real embeddings of K. If this number is positive, then K is formally real and 
consequently its level equals s(K) = oo. Next, Algorithm fl] verifies if —1 is a square 
in K, if so then s(K) = 1. Otherwise, s(K) is either 2 or 4. In order to distinguish 
between these two cases, pj Proposition XI.2.11] comes in handy. It asserts that 
s(K) = 4 if and only if there is 1 < j < k such that dj = ejfj is odd. Otherwise, 
s(K) = 2. This is precisely what step Hat the end of the algorithm is for. □ 

The pseudo-code of this algorithm is presented in Listing [10] 

6. Witt equivalence 

The ultimate problem in the algebraic theory of quadratic forms is to find criteria 
for an existence of an isomorphism between the Witt rings of two fields. Such fields 
are then called Witt equivalent if the above-mentioned isomorphism exists. In this 
section we present an algorithm computing the complete set of Witt equivalence 
invariants of a given number fields. In particular, comparing the results returned 
by the algorithm one can check whether two number fields are Witt equivalent or 
not. It is known (see e.g. [T^]) that the following invariants fully determine the 
Witt class of a number field K: 
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• d-(K:Q) the degree of K over Q; 

• r the number of real embeddings of K: 

• s = s(K) the level of K; 

• k the number of dyadic primes of K\ 

• for each dyadic prime 1)j with 1 < j < k, the pair {dj,Sj) consisting of a 
local degree dj — (K^ )j : Q 2 ) and the local level Sj = s(.K" 0i ). 

We claim that all these invariants are computable. 

Let again K = Q($) be a fixed number field specified by the minimal polyno- 
mial / £ Q[x] of the generator ■&. The first two invariants d and r are trivially 
computable. The degree d is just the degree deg / of the defining polynomial. In 
order to compute r one simply counts the number of real roots of / using Sturm's 
sequence (see e.g. [HI §8.4]). In the previous section we showed how to compute the 
level of K. This leaves us only with the local invariants. Assume that the principal 
ideal 20k factors into prime ideals as: 

20 k = Of • • • d e k k 

and let fj — (Ok/Dj : F 2 ) be the inertia degree of 0^ (1 < j < k). The local degree 
dj = (Kj, : Q 2 ) is the product dj — ejfj. What we need is to determine the local 
level Sj = s(K$.). Fix an even prime = dj. 

Algorithm 13. Let be an even prime of a number field K, e be the ramification 
and / the inertia degree of 0. This algorithm computes the level s(Kj,) of the dyadic 
completion Kx, of K. 

(1) Use Algorithm fl] to check if — 1 is a square in K, if so then return 1 and 
quit. 

(2) If both e and / are odd, then return 4. 

(3) If e is odd but / is even, then return 2. 

(4) If e is even, use Algorithm [3] to check whether —1 is a square in K^,, if so 
then return 1, if not then return 2. 

Proof of correctness. It is clear that if —1 is a square already in K, then it is also a 
square in K a and so s(Ks) — 1. This justifies the first step. Suppose that e is odd. 
Let Q 2 (??) be the (unique) maximal unramified extension of Q 2 contained in Kq. 
Since the quadratic extension Q 2 («)/Q 2 is totally ramified (see e.g. [H Ch. V §2]), 
it follows that i <£ Q 2 (v)- Now, (Q 2 (77) : Q 2 ) = / and (Kj,: Q 2 ) = ef. Hence the 
relative degree (K : Q 2 (?7)) equals e and so is odd. In particular i ^ K- 0l either. 
Thus s(Kq) > 2. Finally Example XI.2.4] asserts that s^K^,) = 4 if and only if 
(K* : Q 2 ) is odd. 

Conversely, assume that e is even and so is the degree (K : Q 2 )- It follows from 
Example XI.2.4] that s(Kx,) < 2. It equals one if and only if —1 is a square 
mi^o. □ 



The pseudo-code of this algorithm is contained in Listing 11 Having all the 
necessary ingredients ready we may now present the last algorithm of this paper 
that construct the complete set of Witt equivalence invariants. 

Algorithm 14. If K = Q(#) is a number field specified (up to an isomorphism) 
by the minimal polynomial / G Q[t] of its generator, then this algorithm computes 
the complete set of Witt equivalence invariants ol K . In particular, two fields are 
Witt equivalent if and only if the outputs of the algorithm are the same for both 
fields. 

(1) Let d = deg/. 

(2) Use Sturm's sequence to compute the number r of real roots of /. 



(3) Use Algorithm 12 to compute the level s = s{K). 
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(4) Use [1, §6.2.5] to factor 20 K , let C = {(Qj, e,-, /_,•)} be the output of this 
algorithm. Take an empty list S. 



(5) For each even prime X>j £ C let dj = &//,-. Use Algorithm 13 to compute 
the local level Sj = s(Kq.). Append the pair [dj,Sj) to the list S. 

(6) Sort the list S lexicographically. 

(7) Return (d,r,s,k,S). 

Observe that the expensive step of decomposing 2 in Ok appears twice: in 
Algorithm [14] as well as in its subroutine computing the level of K. The same can 
be said about reducibility of x 1 + 1 in K[x\. Thus in an actual implementation one 
can combine algorithms 1 1 2[fT4"| into a single procedure compromising the readability. 

7. Example applications 

In order to verify the correctness of the algorithm as well as to allow experi- 
mentation, we implemented the presented algorithm in a computer algebra system 
Sage [11] , A formula for the number of Witt classes of number fields of a fixed 
degree was developed in [12]. Nevertheless, actual representatives for these classes 
was only found for quadratic and cubic fields in 12] and for quartic fields in [6]. 
The first test for usability of our implementation was to find new representatives of 
all 29 classes of quartic fields. It turned out that just by a random search, our im- 
plementation was able to find representatives with generating polynomials having 
smaller coefficients (the biggest coefficient in our case is 752 vs. 208 042 in 6J). The 
results are presented in Table [I] Next, we found the representatives of all 36 classes 
of quintic fields (see Table [2]) and all 95 classes of sextic fields (see Table [3]) . These 
two results are completely new, showing the usability of the developed algorithms. 

Finally, we used the implementation to test empirically how the different Witt 
classes are distributed. To this end, for each of the degrees 3, 4 and 5, we randomly 
generated 5000 fields and counted the numbers of different Witt classes among the 
outputs. Tables HHH] and charts [TH3] summarize the findings. 



Remark 7.1. Observe that some of Witt classes are extremely rare and virtually 
impossible to be found by a blind random search (in particular, such a blind search 



was not able to find a single representatives of classes ( xxxvi ) , ( xxiv ) , ( xii ) of quintic 
fields, even after 24 hours of computing time). These are mostly the classes of fields 
were 2 splits completely. In order to find the representatives of these classes one 
may proceed as follows. Denote 1 0J2 := 2~ ord2a the canonical dyadic norm and 
let || (a-o, ■ • ■ ,»d)||2 '■— max{|ao|2, ■ ■ ■ , |a<i|d} be the associated norm of the vector 
space Ql2) d+l ■ Take a polynomial / with d distinct integral roots. For example, for 
quintic fields we used just 

/ = (x - l)(x - 3) (as - 5)(x - 7) (a; - 11). 

Write it as a dot product f = V ■ X, where V is the vector of coefficients of / and 
X = (1, x, x 2 , . . . , x d ) T are the powers of x. Take now some random vector W and 
let / = (V + W) ■ X. If the norm ||W|| 2 of W is small enough, then / still has 
d distinct roots in Q 2 but there is a good chance taht it is irreducible over Q. It 
follows that 2 splits completely in the field K = Q[x]/(f), as desired. 
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Table 1. Witt classes of quartic fields 



No. 


defining polynomial 


r 


s 


dyadic degrees and levels 


i 


x A + 2x 3 - x 2 - 2x + 2 





1 


{(2,1), (2,1)} 


ii 


x 4 + l 





I 


{(4,1)} 


iii 


x 4 + 2x 3 + 27x 2 + 2x + 2 





2 


{(2,1), (2,1)} 


iv 


x 4 + x 3 + x 2 - 4x + 2 





2 


{(2,1), (2,2)} 


V 


X — X 3 + X + 1 





2 


{(2,2), (2,2)} 


vi 


x 4 - 2x 3 - 3x 2 + 4x + 9 





2 


{(4,1)} 


vii 


x 4 + X + 1 





2 


{(4,2)} 


viii 


x 4 - 29x 3 + 402.x 2 + 752x + 608 





4 


{(1,4),(1,4),(1,4),(1,4)} 


ix 


x 4 - x 3 + f Ox 2 - 2x + 4 





4 


{(1,4), (1,4), (2,1)} 


X 


x 4 + x + 2 





4 


{(1,4), (1,4), (2,2)} 


xi 


x 4 + x 2 — x + 1 





4 


{(1,4), (3,4)} 


xii 


x 4 + 8x 3 ~ 83x 2 - 70a; - 96 


2 


+oo 


{(1,4), (1,4), (1,4), (1,4)} 


xiii 


x 4 - 6x 3 + x 2 - 2x - 2 


2 


+oo 


{(1,4), (1,4), (2,1)} 


xiv 


x 4 - 2x 2 + x - 2 


2 


+oo 


{(1,4), (1,4), (2,2)} 


XV 


x 4 + x 3 - x - 2 


2 


+oo 


{(1,4), (3,4)} 


xvi 


x 4 + 6x 3 - x 2 + 2x + 1 


2 


+oo 


{(2,1), (2,1)} 


xvii 


x 4 - x 3 - x 2 - 2 


2 


+oo 


{(2,1), (2,2)} 


xviii 


x 4 + x 3 + x 2 - 2 


2 


+oo 


{(2,2), (2,2)} 


xix 


,Xj k)JU ZjJu O 


2 


+oo 


{(4,1)} 


XX 


x 4 — x 2 — \ 


2 


+oo 


{(4,2)} 


xxi 


x 4 - 92a; 3 + 39a; 2 + 56a; + 4 


4 


+oo 


{(1,4), (1,4), (1,4), (1,4)} 


xxii 


x 4 - 'ix 3 - f Ox 2 + f 0a; + 4 


4 


+oo 


{(1,4), (1,4), (2,1)} 


xxiii 


x 4 + 2x 3 - f Ox 2 - 5x + 2 


4 


+oo 


{(1,4), (1,4), (2,2)} 


xxiv 


x 4 + x 3 - 4x 2 - x + 2 


4 


+oo 


{(1,4), (3,4)} 


XXV 


x 4 - f Ox 3 - 5x 2 + 6x + 2 


4 


+oo 


{(2,1), (2,1)} 


xxvi 


x 4 — x 3 — 6x 2 — x + f 


4 


+oo 


{(2,1), (2,2)} 


xxvii 


x 4 + 3x 3 - 4x 2 - 3x + f 


4 


+oo 


{(2,2), (2,2)} 


xxviii 


x 4 + 4x 3 - 4x + f 


4 


+oo 


{(4,1)} 


xxix 


x 4 + x 3 - 3x 2 - x + f 


4 


+oo 


{(4,2)} 
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Table 2. Witt classes of quintic fields 



1:5 



No. 


defining polynomial 


r 


dyadic degrees and levels 


i 


X 5 + x 2 + 1 


1 


{(5,4)} 


ii 


x 5 - 2x A + x 3 + x + 2 


1 


{(M), (4,1)} 


iii 


x 5 + x 4 + x 2 + 1 


1 


{(1,4),(4,2)} 


iv 


x 5 — x 3 + x 2 - 2 


1 


{(2,1), (3, 4)} 


V 


x 5 - x + 1 


1 


{(2,2), (3,4)} 


vi 


x 5 - x 3 + x 2 — x + 2 


1 


{(1,4),(1,4),(3,4)} 


vii 


x 5 + 3x 3 - 2x - 4 


1 


{(1,4), (2,1), (2,1)} 


viii 


x 5 - 2a; 4 - x 2 + 2x - 2 


1 


{(1,4), (2,1), (2, 2)} 


ix 


x 5 - x 2 - 2 


1 


{(1,4), (2,2), (2,2)} 


X 


x 5 - 4x 4 - 3x 3 + 2x - 4 


1 


{(1,4), (1,4), (1,4), (2,1)} 


xi 


x 5 - x 4 - 2x 3 - 5x 2 + x - 2 


1 


{(1,4), (1,4), (1,4), (2,2)} 


xii 


x 5 + 37x 4 + 774x 3 - f 058x 2 + 2f 69x - 643 


1 


{(1,4),(1,4),(1,4),(1,4),(1,4)} 


xiii 


X 5 + x 4 + x 2 — x — 1 


3 


{(5,4)} 


xiv 


x 5 + x 4 - x 3 — x 2 - 3.T + 1 


3 




'(1.4), (4,1)} 


XV 


x 5 - 2x 3 - 2a; - 1 


3 




(1,4), (4,2)} 


xvi 


x 5 - x 4 - 2x 3 - x 2 - 2x + 2 


3 




(2,1), (3,4)} 


xvii 


x 5 - x 4 - 2x 3 + 1 


3 


{(2.2).(3.4)} 


xviii 


x 5 - x 3 + x 2 - 3x - 2 


3 




;(1,4),(1,4),(3,4)} 


xix 


x 5 - bx 3 - 2x - 4 


3 




^1,4), (2,1), (2,1)} 


XX 


x 5 - 2x 4 - 3x 3 + 2x + 4 


3 




(1,4), (2,1), (2, 2)} 


xxi 


x 5 + x 4 + x 2 - 'ix - 2 


3 


{(1,4), (2, 2), (2, 2)} 


xxii 


x 5 + 26x 4 + f 9x 3 + 24x 2 - 26x - f 2 


3 


{(1,4), (1,4), (1,4), (2,1)} 


xxiii 


x 5 + fOx 4 - fix 2 - 34x - 16 


3 


{(1,4), (1,4), (1,4), (2,2)} 


xxiv 


x 5 + f Ofx 4 + 326x 3 + 926x 2 + 2f 05x + 893 


3 


{(1,4),(1,4),(1,4),(1,4),(1,4)} 


XXV 


x 5 + 4x 4 - 5x 2 + f 


5 


{(5,4)} 


xxvi 


x 5 + x 4 - f fx 3 - x 2 + 7x + f 


5 




'(1,4), (4,1)} 


xxvii 


x 5 - 5x 4 + f 3x 2 - 5 


5 




> (1,4),(4,2)} 


xxviii 


x 5 - f 5x 4 + f 8x 3 + 37x 2 - f 6x - f 2 


5 




(2,1), (3, 4)} 


xxix 


x 5 - 7x 4 + 6x 2 - f 


5 




(2,2), (3,4)} 


XXX 


x 5 + x 4 - 9x 3 - f 2x 2 + x + 2 


5 




'(1,4), (1,4), (3,4)} 


xxxi 


x 5 + 3f x 4 + 244x 3 + 524x 2 + 99x - f 9 


5 




> (1,4),(2,1),(2,1)} 


xxxii 


x 5 + 26x 4 - 5x 3 - 93x 2 + 42x + f 7 


5 




(1,4), (2,1), (2, 2)} 


xxxiii 


x 5 +x 4 - ffx 3 + f4x-4 


5 




(1,4), (2,2), (2,2)} 


xxxiv 


x 5 + 247x 4 + 2064x 3 + 41 44x 2 + 1887x + 233 


5 


{(1,4), (1,4), (1,4), (2,1)} 


XXXV 


x 5 - 9x 4 + f 4x 3 + 29x 2 - 21 x + 2 


5 


{(1,4), (1,4), (1,4), (2,2)} 


xxxvi 


x 5 + 997x 4 + 33030x 3 + 31646x 2 + 6137x - 131 


5 


{(1,4),(1,4),(1,4),(1,4),(1,4)} 
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Table 3: Witt classes of sextic fields 



No. 


defining polynomial 


r 


s 


dyadic degrees and levels 


i 


x b + x 4 - 2x 3 + 1 





1 


{(6,1)} 




ii 


x 6 + x 4 + 2x 3 + 2x 2 + 2 





I 


{(2,1), (4,1)} 




iii 


x 6 +5x 4 -f6x 3 +4x 2 +32x+64 





I 


{(2,1), (2,1), (2,1)} 




iv 


x 6 + x 4 + 2x 3 + 8x 2 + f 





2 


{(6,1)1 




V 


x 6 + x + f 





2 


1(6,2)} 




vi 


x 6 + x 4 + 15x 2 + 2x + 2 





2 




;(2,1),(4,1)} 




vii 


x 6 + x 3 - x 2 + 2 





2 




(2,1), (4,2)} 




viii 


x 6 + 2x 5 + x 4 - 2x 3 + 2x 2 + 2 





2 




(2, 2), (4,1)} 




ix 


x 6 + x 3 + x 2 + x + f 





2 


{(2,2),(4,2)} 




X 


x 6 + 15x 4 - 30x 3 + 86x 2 - 
72x + 72 





2 


{(2,1), (2,1), (2,1)} 




xi 


x 6 - x 5 + 6x 4 - x 3 - x 2 + 2 





2 




[(2, 1), (2,1), (2,2)} 




xii 


x 6 - x 5 + x 3 - x 2 + 2 





2 




(2,1), (2, 2), (2, 2)} 




xiii 


x 6 + x 5 - x 3 + x 2 + 2 





2 




(2, 2), (2, 2), (2, 2)} 




xiv 


x 6 + x 3 + x 2 + f 





4 


{(1,4), (5,4)} 




XV 


x 6 + x 5 - x 3 + f 





4 


{(3,4), (3,4)} 




xvi 


x 6 - x 5 - x 4 + x 3 + x 2 — x + 2 





4 


{(1,4),(1,4),(4,1)} 




xvii 


x 6 - x + 2 





4 


{(1,4),(1,4),(4,2)} 




xviii 


x 6 + x 4 - x 3 + x 2 + 2 





4 


{(1,4), (2,1), (3,4)} 




xix 


x 6 + x 5 + x 2 + f 





4 


{(1,4), (2, 2), (3,4)} 




XX 


x 6 - x 5 + 5x 4 - 2x 2 + x + 4 





4 


{(1,4), (1,4), (1,4), (3,4)} 




xxi 


x 6 + 8x 5 + 26x 4 + 30x 3 + 
33x 2 + 6x + 8 





4 


{(1,4), (1,4), (2,1), (2,1)} 




xxii 


x 6 + x 5 - 3x 3 + 17x 2 + 4 





4 


{(1,4),(1,4),(2,1),(2,2)} 




xxiii 


x 6 + x 4 + 3x 3 + 2x 2 - x + 2 





4 


{(1,4),(1,4),(2,2),(2,2)} 




xxiv 


x 6 + 36x 5 + f Oflx 4 + 
224x 3 + 95x 2 - 4x + 173 





4 


{(1,4), (1,4), (1,4), (1,4), (2,1)} 




XXV 


x 6 + If 2x 5 + 8265x 4 + 160x 3 + 
7179x 2 + 2544x + 273f 





4 


{(1,4), (1,4), (1,4), (1,4), (2,2)} 




xxvi 


x 6 + f 9x 4 - f 6x 3 + f48x 2 + 
352x + 512 





4 


{(1,4), (1,4), (1,4), (1,4), (1,4), (1,4 


)} 


xxvii 


x 6 + x 4 - 2x 3 - 3 


2 


+oo 


{(6,1)} 




xxviii 


x 6 + x 5 - 1 


2 


+oo 


{(6,2)} 




xxix 


x 6 — x 5 — x 3 — x 2 — X + f 


2 


+oo 


{(1,4), (5,4)} 




XXX 


x 6 - 4x 5 - 2x 4 + 2x 3 - 2x 2 - 
2x + f 


2 


+oo 


{(2,1), (4,1)} 




xxxi 


x 6 + x 5 + x 4 + x 3 - x 2 - 2 


2 


+oo 




'(2,1), (4, 2)} 




xxxii 


x 6 -3 


2 


+oo 




(2,2), (4,1)} 




xxxiii 


x 6 — x 5 — x — f 


2 


+oo 




,(2,2), (4,2)} 




xxxiv 


X 6 + X 5 + x 3 — f 


2 


+oo 




(3,4), (3,4)} 




XXXV 


x 6 - x 5 - x 4 + x 3 + x 2 + x - 4 


2 


+oo 


\ 


;(1,4),(1,4),(4,1)] 






xxxvi 


x 6 — x 4 — x 2 — x — 2 


2 


+oo 




'(1,4), (1,4), (4,2)' 






xxxvii 


x 6 + 2x 5 - x 4 - x 3 - x 2 - 2 


2 


+oo 




^1,4), (2,1), (3,4)' 






xxxviii 


x 6 - x 3 + x 2 — x — 2 


2 


+oo 




(1,4), (2, 2), (3, 4)' 






xxxix 


x 6 - 2x 5 - 2x 4 + f Ox 3 + 


2 


+oo 




(2,1), (2,1), (2,1) 








25x 2 - 28x - 28 










xl 


x 6 + 3x 5 - x 3 + x 2 - 2 


2 


+oo 


{(2,1), (2,1), (2,2)} 




xli 


x 6 - x 5 + x 3 + x 2 - 2x - 2 


2 


+oo 


{(2,1), (2, 2), (2, 2)} 





Continued on the next page 
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Table 3: Witt classes of sextic fields 



No. 


defining polynomial 


r 


s 


dyadic degrees and levels 


xlii 


x B -x 5 -2x 4 -x 3 + x 2 -2x-2 


2 


+oo 


{(2, 2), (2, 2), (2, 2)} 


xliii 


x 6 + 5x 5 + x 4 + 2x 3 + x - 2 


2 


+oo 


{(1,4), (1,4), (1,4), (3,4)} 


xliv 


x 6 - 6x 4 - 2x 3 + x 2 - 2x - 8 


2 


+oo 


{(1,4), (1,4), (2,1), (2,1)} 


xlv 


x 6 - 2x 5 - 2x 4 + x 3 - 2x 2 + 
2x -8 


2 


+oo 


{(1,4), (1,4), (2,1), (2,2)} 


xlvi 


x 6 -3x 5 -x 4 + x 3 -x 2 +x-2 


2 


+oo 


{(1,4), (1,4), (2,2), (2,2)} 


xlvii 


x 6 - 232a; 5 - 479x 4 - 
440a; 3 - 502a; 2 + 348x - 64 


2 


+oo 


{(1,4), (1,4), (1,4), (1,4), (2,1)} 


xlviii 


x 6 - 45x 5 - 38a; 4 - x 3 + 
41a; 2 - 6a; - 16 


2 


+oo 


{(1,4), (1,4), (1,4), (1,4), (2,2)} 


xlix 


x 6 - 53a; 4 - 16a; 3 + 868a; 2 - 
800a; - 4288 


2 


+oo 


{(1,4), (1,4), (1,4), (1,4), (1,4), (1,4)} 


1 


x 6 + 6a; 5 + 2a; 4 - 2a; 3 - 5x 2 - 
2a; + 1 


4 


+oo 


{(6,1)} 


li 


x 6 - x 5 - 2x 4 - x 2 + x + 1 


4 


+oo 


{(6,2)} 


lii 


x 6 + x 5 - 2x 4 - x 3 - x 2 - x + 2 


4 


+oo 


\ 


'(1,4),(5,4)} 


liii 


x 6 + 8x 5 - a; 4 - 4a; 3 - 4a; 2 + 2 


4 


+oo 




(2,1), (4,1)} 


liv 


x 6 - a; 5 - x 4 + x 3 - 3a; 2 + 2 


4 


+oo 




'(2,1), (4, 2)] 




lv 


x 6 - 3a; 5 +x 4 -3x 2 +x + l 


4 


+oo 




(2,2), (4,1)' 




lvi 


x 6 -2x 5 -x 4 -x 3 -x 2 + 2x + l 


4 


+oo 




(2, 2), (4, 2) 




lvii 


x 6 + x 5 -x 4 + x 3 - x 2 - 3x + 1 


4 


+oo 




"(3,4), (3, 4) 




lviii 


x 6 - 14a; 5 - 6a; 4 + 4.T 3 - 
6a; 2 + 26.T + 19 


4 


+oo 


{(1,4),(1,4),(4,1)} 


lix 


x 6 + 2a; 5 - 2a; 4 - 2a; 3 - 3a; + 2 


4 


+oo 




•(1,4), (1,4), (4,2)} 


lx 


.x 6 +a; 5 -3x 4 -2a; 3 -2.a; 2 +x+6 


4 


+oo 




(1,4), (2,1), (3,4)} 


lxi 


x e + x 5 - 2x 4 - x 2 - Ax + 1 


4 


+oo 




(1,4), (2, 2), (3, 4)} 


lxii 


x 6 - 2x 5 - 10a; 4 - 2a; 3 - 


4 


+oo 




'(2,1), (2,1), (2,1)} 




llx 2 -16a; + 16 








lxiii 


x 6 + 13a; 5 + 24a; 4 - 47a; 3 - 
37x 2 - 44a; - 46 


4 


+oo 


{(2,1), (2,1), (2,2)} 


lxiv 


.x 6 - 3a; 5 - 4a; 4 + 3a; 3 - x 2 + 2 


4 


+oo 


{(2,1), (2, 2), (2,2)} 


lxv 


x 6 - x 5 - 4x 4 + x 3 - 3x 2 + 2 


4 


+oo 


{(2,2), (2,2), (2,2)} 


lxvi 


x 6 + x 5 + x 4 - 21x 2 + 4 


4 


+oo 


{(1,4), (1,4), (1,4), (3,4)} 


lxvii 


x 6 - 484a; 5 - 696x 4 + 
346a; 3 - 671x 2 - 538x + 114 


4 


+oo 


{(1,4), (1,4), (2,1), (2,1)} 


lxviii 


x 6 - 15x 5 + llx 4 + 25x 3 + 
32x 2 + 34x + 8 


4 


+oo 


{(1,4), (1,4), (2,1), (2,2)} 


lxix 


x 6 + x 5 -4x 4 -x 3 -x 2 -2x + 8 


4 


+oo 


{(1,4), (1,4), (2,2), (2,2)} 


lxx 


x 6 + 474a; 5 - 72x 4 - 424x 3 + 
415x 2 - 434x - 280 


4 


+oo 


{(1,4), (1,4), (1,4), (1,4), (2,1)} 


lxxi 


x 6 + 120a; 5 + 2x 4 - 169x 3 - 
170x 2 + 176a; + 32 


4 


+oo 


{(1,4), (1,4), (1,4), (1,4), (2,2)} 


lxxii 


x 6 + 8152x 5 + 131685x 4 + 
3664x 3 + 49395x 2 + 
496600x + 23207 


4 


+oo 


{(1,4), (1,4), (1,4), (1,4), (1,4), (1,4)} 


lxxiii 


x 6 - 27x 4 - 2x 3 + 156x 2 - 
36x- 111 


6 


+oo 


{(6,1)} 


lxxiv 


x 6 - 10x 4 + x 3 + 8x 2 - 1 


6 


+oo 


{(6,2)} 


lxxv 


x 6 -2x 5 -46x 4 +71x 2 +27x+l 


6 


+oo 


{(1,4), (5,4)} 



Continued on the next page 
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Table 3: Witt classes of sextic fields 



No. 


defining polynomial 


r 


s 


dyadic degrees and levels 


lxxvi 


x B - 29x 4 - 2x 3 + 163x 2 - 
34x - 62 


6 


+oo 


{(2,1), (4,1)} 


lxxvii 


x 6 - 57x 5 - 137x 4 + 325a; 3 + 
557x 2 - 372x - 86 


6 


+oo 


{(2,1), (4, 2)} 


lxxviii 


x 6 -20x 4 -4x 3 +61x 2 +4x-46 


6 


+oo 


{(2,2), (4,1)} 


lxxix 


x 6 - 16x 4 - x 3 + 15x 2 - 3a; - 1 


6 


+oo 


{(2,2), (4,2)} 


lxxx 


x 6 - 42a; 5 + f fa; 4 + 414a; 3 - 
464a; 2 - 386x + 423 


6 


+oo 


{(3,4), (3,4)} 


lxxxi 


x e - 270a; 5 - 478a; 4 + 
836a; 3 + 890x 2 - 252x + 9 


6 


+oo 


{(1,4), (1,4), (4,1)} 


lxxxii 


x 6 - 69a; 5 + 320a; 4 + 388a; 3 - 
252a; 2 - 122x + 8 


6 


+oo 


{(1,4), (1,4), (4,2)} 


lxxxiii 


x 6 - f 40a; 5 + 285a; 4 + 
507a; 3 - 508x 2 - 470x - 64 


6 


+oo 


{(1,4), (2,1), (3,4)} 


lxxxiv 


x e - 2x 5 - 35x 4 + 20a; 2 +x-l 


6 


+oo 


{(1,4), (2,2), (3,4)} 


lxxxv 


x 6 - 35x 4 - 4a; 3 + f 96a; 2 - 
56a; + 4 


6 


+oo 


{(2,1), (2,1), (2,1)} 


lxxxvi 


x 6 + 64a; 5 + f 32a; 4 + 2a; 3 - 
43x 2 + 2x + 2 


6 


+oo 


{(2,1), (2,1), (2,2)} 


lxxxvii 


x 6 + 1 28a; 5 + 2049x 4 + 

81 92a; 3 + 8151x 2 + 128a; - 41 


6 


+oo 


{(2,1), (2,2), (2,2)} 


lxxxviii 


x 6 - 19a; 5 - 34x 4 + 127a; 3 + 
169x 2 - 88x - 54 


6 


+oo 


{(2,2), (2,2), (2,2)} 


lxxxix 


x 6 + 8a; 5 - 105a; 4 + 8x 3 + 
912a; 2 + 16x - 768 


6 


+oo 


{(1,4),(1,4),(1,4),(3,4)} 


xc 


x 6 + 2x 5 - 77x 4 + 4a; 3 + 
910x 2 + 4.x - 768 


6 


+oo 


{(1,4), (1,4), (2,1), (2,1)} 


xci 


x 6 + 4x 5 - 105x 4 + 4x 3 + 
910x 2 + 4x - 768 


6 


+oo 


{(1,4),(1,4),(2,1),(2,2)} 


xcii 


x 6 - 283x 5 + 131x 4 + 
448x 3 - 238x 2 - 28x + 16 


6 


+oo 


{(1,4),(1,4),(2,2),(2,2)} 


xciii 


x 6 + 1028x 5 + 131059x 4 + 
4194272x 3 + 8388639x 2 + 
4194268x + 131117 


6 


+oo 


{(1,4), (1,4), (1,4), (1,4), (2,1)} 


xciv 


x 6 +496x 5 + 65611x 4 + 
2097152x 3 + 2096311x 2 + 
6032x - 643 


6 


+oo 


{(1,4), (1,4), (1,4), (1,4), (2,2)} 


xcv 


x 6 - 85x 4 - 16x 3 + 1156x 2 - 
544x + 64 


6 


+oo 


{(1,4),(1,4),(1,4),(1,4),(1,4),(1,4)} 
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Table 4. Number of different Witt classes in 5000 random cubic fields. 



No. 


r 


dyadic degrees and levels 


# in 5000 
random fields 


i 


1 


{(3,4)} 


1647 


ii 


1 


{(1,4), (2,1)} 
{(1,4), (2,2)} 


71 


iii 


1 


908 


iv 


1 


{(1,4), (1,4), (1,4)} 


23 


V 


3 


{(3,4)} 


470 


vi 


3 


{(1,4),(2,1)} 
{(1,4), (2,2)} 


36 


vii 


3 


202 


viii 


3 


{(1,4), (1,4), (1,4)} 


2 



Table 5. Number of different Witt classes in 5000 random quartic fields. 



No. 


r 


s 


dyadic degrees and levels 


# in 5000 
random fields 


i 





1 


r {(4 ' 1)j , 


46 


ii 





1 


{(2,1), (2,1)} 





iii 





2 


{(4,1)} 


5 


iv 





2 


{(4,2)} 


321 


V 





2 




(2,1), (2,1)} 





vi 





2 




(2,1), (2, 2)} 


16 


vii 





2 




.(2,2), (2,2)} 


108 


viii 





4 




'(1,4), (3,4)} 


315 


ix 





4 


{(1,4), (1,4), (2,1)} 


3 


X 





4 


{(1,4), (1,4), (2,2)} 


48 


xi 





4 


{(1,4), (1,4), (1,4), (1,4)} 





xii 


2 


+oo 


{(4,1)} 


33 


xiii 


2 


+oo 


{(4,2)} 


1087 


xiv 


2 


+oo 




(1,4), (3,4)} 


804 


XV 


2 


+oo 




(2,1), (2,1)} 


2 


xvi 


2 


+oo 




(2,1), (2, 2)} 


85 


xvii 


2 


+oo 


{ 


'(2,2), (2,2)} 


287 


xviii 


2 


+oo 


{(1,4), (1,4), (2,1)} 


12 


xix 


2 


+oo 


{(1,4), (1,4), (2,2)} 


121 


XX 


2 


+oo 


{(1,4), (1,4), (1,4), (1,4)} 





xxi 


4 


+oo 


{(4,1)1 


6 


xxii 


4 


+oo 


{(4,2)} 


82 


xxiii 


4 


+oo 




'(1,4), (3,4)} 


48 


xxiv 


4 


+oo 




(2,1), (2,1)} 





XXV 


4 


+oo 




(2,1), (2, 2)} 


2 


xxvi 


4 


+oo 




'(2,2), (2,2)} 


16 


xxvii 


4 


+oo 


{(1,4), (1,4), (2,1)} 





xxviii 


4 


+oo 


{(1,4), (1,4), (2,2)} 


3 


xxix 


4 


+oo 


{(1,4), (1,4), (1,4), (1,4)} 
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Table 6. Number of different Witt classes in 5000 random quintic fields. 



No. 


r 


i 


f 


ii 


f 


iii 


f 


iv 


f 


V 


f 


vi 


f 


vii 


f 


viii 


f 


ix 


f 


X 


f 


xi 


f 


xii 


f 


xiii 


3 


xiv 


3 


XV 


3 


xvi 


3 


xvii 


3 


xviii 


3 


xix 


3 


XX 


3 


xxi 


3 


xxii 


3 


xxiii 


3 


xxiv 


3 


XXV 


5 


xxvi 


5 


xxvii 


5 


xxviii 


5 


xxix 


5 


XXX 


5 


xxxi 


5 


xxxii 


5 


xxxiii 


5 


xxxiv 


5 


XXXV 


5 


xxxvi 


5 



dyadic degrees and levels 



# in 5000 
random fields 



(I- 

(1, 

"(1, 

{(1. 
(1,4), 

'(1,4), 

{(M),(i, 



(1, 
(i, 
(i, 
(i, 

(1,4), 
'(1,4), 

{(1,4), (f, 



{(1, 
"(1, 
(1, 
(1, 
{(1,4), 
{(1,4), 
{(1,4), (1, 



{(5,4)} 

(1,4), (4, f) 
(1,4), (4, 2) 
(2,1), (3,4) 
(2, 2), (3, 4) 
4), (1,4), (3 
4), (2,1), (2 
4), (2,1), (2 
4), (2, 2), (2 
(1,4), (1,4) 
(1,4), (1,4) 
4),(1,4),(1 

{(5,4)} 
(1,4), (4,1) 
(1,4), (4, 2) 
(2,1), (3,4) 
(2, 2), (3, 4) 
4), (1,4), (3 
4), (2,1), (2 
4), (2,1), (2 
4), (2, 2), (2 
(1,4), (1,4) 
(1,4), (1,4) 
4),(1,4),(1 

{(5,4)} 
(1,4), (4,1) 
(1,4), (4, 2) 
(2,1), (3,4) 
(2, 2), (3, 4) 
4), (1,4), (3 
4), (2,1), (2 
4), (2,1), (2 
4), (2, 2), (2 
(1,4), (1,4) 
(1,4), (1,4) 
4),(1,4),(1 




4) 

1) 

2) 

2)' 

(2,1) 

(2,2) 

4), (1,4)} 




768 

25 

598 

53 

513 

132 



24 

91 



7 



497 

17 

311 

33 

292 

53 

3 

10 

43 



2 



6 

1 

4 



5 

















COMPUTING WITH QUADRATIC FORMS OVER NUMBER FIELDS 



J. i) 



1600 



1200 



800 



400 





12 3 4 5 6 7 



Figure 1. Number of different Witt classes in 5000 random cubic fields. 



1000 



300 



600 



400 



200 



I 






1 1 



1 



7 



10 13 16 19 22 25 



28 



Figure 2. Number of different Witt classes in 5000 random quar- 
tic fields. 
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Figure 3. Number of different Witt classes in 5000 random quin- 
tic fields. 
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Appendix A. Pseudo-codes of presented algorithms 

For convenience of the readers, especially those wishing to implement the pre- 
sented procedures, this section gathers the pseudo-codes of all described algorithms. 

Input: a and element of a number field K 

{true if a is a square in K 
false otherwise 

// Use [15, §5.4] 

if is_irreducble(a; 2 — a) then 
return true; 

return false; 

Listing 1: is_square 



Input: < 
Output: 



a a non-zero element K 

p an odd prime of a number field K 

{true if a is a square in K p 
false otherwise 



if ordp a = 1 (mod 2) then 
return false; 

// Use Jl], Algorithm 4.8.17] 

7r <— unif ormizer(p); 

a' <r- a-ir- m ' d r a ; 

if is_square(a', K/p) then 

return true; 
else 
_ return false; 

Listing 2: is_nondyadic_square 



Input: ■ 
Output 



a a non-zero element K 

d an even prime of a number field K 

{true if a is a square in K^, 
false otherwise 

// Use Algorithm fl] 

if is_square(a; K) then 
return true; 

L^- K[x]/(x 2 -a); 

// Use O, Algorithm 2.4.13] 

C <— prime_decomposition(c); L); 

if jl£ = 2 then 
return true; 
else 
[_ return false; 

Listing 3: is_dyadic_squarc 



Output: 
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{tp = (oi, . . . , ad) a non-degenerate quadratic form over K, 
all at are p-adic units 
p an odd prime of a number field K 

{true if 93 ® K/p is isotropic 
false otherwise 

switch dim((p) do 
case 1 
I return false; 

case 2 

if is_square(disc((,c), K/p) then 

I return true; 
else 

|_ return false; 

otherwise return true; 

Listing 4: is_residually .isotropic 

{ip = (01, . . . , ad) a non-degenerate quadratic form over K 
p an odd prime of a number field K 

{true if ip ® if p is isotropic 
false otherwise 

if dim(^) = 1 then 
_ return false; 

if dim(< / o) > 5 then 
return true; 

// Use EL Algorithm 4.8.17] 

7r <— unif ormizer(p); 

ipo <- {a, ; • 7r~ 

ipi 4- {c 

II Use Algorithm H 

^4 <— is_residually_isotropic(y>o; p); 

_B <— is_residually_isotropic(y?i;p); 

return 4VB; 

Listing 5: is_locally_isotropic_non-dyadic 



: a, ■ 7r- ord p 



ordp Oj = (mod 2)}; 
ordp Oj = 1 (mod 2) j; 
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= (oi, . . . , ad) a non-degenerate quadratic form over K 
an even prime of a number field K 

{true if ip (g) if is isotropic 
false otherwise 

switch dim(tp) do 
case 1 
|_ return false; 

case 2 

// Use Algorithm [3] 

if is_dyadic_square(— aia 2 , 0) then 

I return true; 
else 

I return false; 

case 3 

// Use Eq. Q and El, Algorithm 6.6] 
S <— rXt< ,- hilbert_symbol(aj, a^; 0); 
// Use [0, Proposition V.3.22] 

if s = hilbert_symbol(— 1, — det(< / 9); c)) then 

I return true; 
else 

I return false; 

case 4 

// Use Algorithm [3] 

if is_dyadic_square(det t/j,i)) then 

// Use Eq. and EH, Algorithm 6.6] 

s <— Yii< j hilbert_symbol(ai, aj;T>); 

II Use [7, Proposition V.3.23] 

if s = hilbert_symbol(— 1, —1; d) then 

I return true; 
else 

|_ return false; 

else 
I return true; 

otherwise return true; 

Listing 6: isdocally_isotropic_dyadic 
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Input: ip = (ai, . . . , a,d) a non-degenerate quadratic form over K 

{true if ip is isotropic 
false otherwise 

if dim(< / o) = 1 then 
_ return false; 

if dim(y) = 2 then 

// Use Algorithm fl] 

if is_square(— axa^; 2<Q then 
|_ return true; 

return false; 

// Check all real places 
C <— real.embeddings^); 
for p e C do 

S <- sgnp(^); 

if \s\ = dimip then 
|_ return false; 

// Check all dyadic primes 
C <— prime_decomposition(2, K); 
for t) € £ do 

// Use Algorithm [6] 

if not is_locally_isotropic_dyadic((/3; d) then 
\_ return false; 

// Find all odd primes that matter 

for 1 < i < dim ip do 

C <— C U prime_decomposition(ai; K); 

C •<— C, \ prime_decomposition(2, K)] 
II Check the odd primes 
for p e C do 

// Use Algorithm [5] 

if not is_locally_isotropic_non-dyadic(<^; p) then 
I return false; 

return true; 

Listing 7: is_globally .isotropic 
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Input: ip = (ai, . . . , a<j) a non-degenerate quadratic form over K 

{true if ip is hyperbolic 
false otherwise 

if dim ip = 1 (mod 2) then 
_ return false; 

// Use Algorithm fl] 

if not is_square(disc <p\ K) then 
return false; 

// Check all real places 
C 4— real.embeddings^); 
for p e C do 

S <- sgnp(^); 

if s 7^ then 
|_ return false; 

// Find all odd primes that matter 

for 1 < i < dimip do 

C <— C U prime_decomposition(ai; K)\ 

C <— C \ prime_decomposition(2, K); 
II Check the odd primes 
for p e C do 

// Use [1, Algorithm 4.8.17] 

7r <— unif ormizer(p); 

(^i <- {a 4 • 7r-° rd > a < | ordp a t = 1 (mod 2)}; 

if dim^i = 1 (mod 2) or not is_square(discipi; K/p) then 
| return false; 

// Check all dyadic primes 
C <— prime_decomposition(2, K); 
for 3e£ do 

// Use Eq. ([!]) and [Q3L Algorithm 6.6] 

s <— J| i<: hilbert_symbol(ai, a,-; 0); 

if s ^ 1 then 
|_ return false; 

return true; 

Listing 8: is_globally .hyperbolic 
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Input: tp = (ai, . . . , ad) a non-degenerate quadratic form over K 
Output: d the dimension of the anisotropic part of p 

II Find the signatures 
C 4— real_embeddings(/ - iT); 
N 4— max{sgn p(p) | p G £}; 
if TV > 3 then 
_ return N; 

ifiV = 0V7V = 2 then 
if N = then 

// Use Algorithm [SI 

if is_globally_hyperbolic(</?) then 
|_ return 0; 

// Find all primes that matter 

for 1 < i < dim <p do 

| C 4— C U prime_decomposition(ai; K); 

C 4— C U prime_decomposition(2, if); 
// Check these primes 
for p e £ do 

// Use Eq. Q and 03], Algorithm 6.6] 

s ■<— Oi< ,- hilbert_symbol(aj, a_j-;p); 

if s 7^ 1 then 
|_ return 4; 

if iV = 1 then 

ip 4— ip -L (det </?) ; 

return (dimension_of _anisotropic_part(?/;) — 1); 

Listing 9: dimension_of_anisotropic_part 

Input: / £ Q[i] the minimal polynomial of the generator ■& 
Output: s = s(K) the level of the number field K = Q(i?) 

// Use Sturm sequence to compute the number of real roots of /: 
r 4— count_real_roots(/); 
if r^O then 
_ return oo; 

// Use Algorithm fl] 

if is_square( — 1; K) then 
return 1; 

// Use [I], §6.2.5] to decompose 2 in Ok, 
II C consists of triples (ftj,ei,/j). 
C 4— prime_decomposition(2, K); 
II Use IH, Proposition XI. 2. 11] 
for (D, e, /) G £ do 

if (e • /) = 1 (mod 2) then 
[_ return 4; 

return 2; 

Listing 10: level 
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{0 an even prime of k 
e the ramification index of p 
/ the inertia degree of p 
Output: s = s(if ) level of the completion of K at D 

II Use Algorithm fl] 

if not is_square(— 1, K) then 
[_ return 1; 

if e = 1 (mod 2) then 
if / = 1 (mod 2) then 

| return 4; 
else 

[_ return 2; 

else 

// Use Algorithm [5] 

if is_dyadic_square(— 1, D) then 

| return 1; 
else 

|_ return 2; 

Listing 11: dyadic Jevel 

Input: / 6 Q[i] the minimal polynomial of the generator of a number field K 
Output: The complete set of Witt equivalence invariants of K in the form 
(d, r, s, k, {(di, si), . . . .(d k , s k )}) 

II The degree is just (K: Q) = deg/. 

d^degf; 

II Use Sturm sequence to compute the number of real roots of /. 

r <— count_real_roots(/); 



// Use Algorithm 12 to compute the level of K. 

s <r- level(/); 

// Use IfT], §6.2.5] to decompose 2 in the maximal order of K. 

II C consists of triples (pi,ej,/j) 

C <— prime_decomposition(2, K); 

II The number of dyadic primes: 

k <- ft£; 

for 1 < i < k do 

// The local degree: 



II Use Algorithm 13 to compute the local level. 

Si <- dyadicaevel^p^e^/i); 

5<-5U{(4«i)}i 
// Sort the dyadic degrees and levels lexicographically, 
sort (5); 
return (d, r, s, k 7 S); 

Listing 12: Witt_equivalence_invariants 



